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ABSTRACT 


The  theory  of  block  copolymer  micelle  formation  formulated 
by  Leibler,  Orland  and  Wheeler  is  compared  with  the  experimental 
results  of  Rigby  and  Roe  obtained  with  mixtures  of 
styrene-butadiene  diblock  copolymer  and  butadiene  homopolymer. 

To  facilitate  the  comparison,  the  theoretical  expressions  are 
recast  in  terms  of  the  molecular  volume  v  and  the  rms  end-to-end 
distance  r,  thereby  avoiding  the  use  of  the  concept  of  the 
segment  size  a  and  the  number  of  segment  N  per  molecule.  The 
calculated  radii  of  the  micelles  show  excellent  agreement  with 
the  observed  values.  The  theory  predicts  correctly  the  trend  of 
change  in  the  critical  micelle  concentration  with  the  change  in 
temperature  and  in  the  relative  lengths  of  the  styrene  and 
butadiene  blocks.  The  theory  predicts  the  micelle  cores  to 
consist  predominantly  of  styrene  blocks,  whereas  the  experimental 
results  indicate  that  above  a  certain  temperature  the  micelle 
cores  become  increasingly  swollen  with  poly  butadiene  as  the 
temperature  is  raised  toward  the  dissolution  temperature. 


INTRODUCTION 

1  P 

In  the  previous  two  papers,  ’  hereafter  referred  to  as 
Part  I  and  Part  II,  we  presented  the  result  of  investigation,  by 
small-angle  X-ray  scattering,  of  the  characteristics  of  spherical 
micelles  formed  when  small  amounts  of  styrene-butadiene  block 
copolymer  is  dissolved  in  poly  butadiene.  The  parameters 
evaluated  include  the  critical  micelle  concentration  (CMC),  the 
radius  of  the  micelle  core,  and  the  degree  of  swelling  of  the 
latter  as  a  function  of  the  temperature,  the  concentration  of  the 
block  copolymer  and  the  relative  lengths  of  the  two  blocks  in  the 
copolymer . 

Theoretical  investigation  of  block  copolymer  micelle 
formation  in  the  medium  of  homopolymer  was  recently  carried  out 
by  Leibler,  Orland  and  Wheeler^  and  by  Noolandi  and  Hong.1*  The 
purpose  of  the  present  paper  is  to  compute  numerical  results 
based  on  the  theory  of  Leibler  et  al.^  and  compare  them  with  the 
experimental  results  presented  in  Parts  I  and  II. 

THEORETICAL  EXPRESSIONS 

Before  setting  out  to  make  numerical  computations,  we  at 
first  introduce  some  minor  modifications  to  the  original 
expressions  given  in  the  paper^  by  Leibler,  Orland  and  Wheeler 
(hereafter  referred  to  as  LOW).  The  modifications  are  of  two 
kinds.  The  first  is  introduced  in  order  to  avoid  the  use  of  the 
concept  of  a  "segment,"  which  is  often  defined  ambiguously.  The 
second  is  to  allow  slight  modification  of  the  physical  model  on 
which  the  theory  is  based  so  that  the  possibility  of  swelling  of 
the  micelle  core  can  be  investigated. 


In  the  theories  in  which  configurational  properties  of 
polymer  chains  enter  into  consideration,  a  polymer  molecule  is 
often  represented  by  a  chain  of  N  segments,  each  of  length  a. 
The  exact  size  of  the  length  parameter  a  is  usually  undefined, 
and  this  poses  no  problem  as  far  as  theoretical  development  is 
concerned.  In  fact,  it  can  be  shown  in  a  general  term  that,  when 
the  theory  deals  only  with  bulk  properties  and  no  surface 
phenomena  are  involved,  the  segment  a  can  be  chosen  entirely 
arbitrarily  without  affecting  the  numerical  predictions  of 
thermodynamic  functions  by  the  theory.  When  the  theory  includes 
the  consideration  of  surface  phenomena,  however,  the  numerical 
results  obtained  depend  on  the  precise  definition  of  a,  as  the 
depth  of  the  surface  effect  predicted  is  proportional  to  the 
segment  size  a.  The  customary  practice  is  to  equate  a  to  some 
dimension  (e.g.,  length  or  cube  root  of  volume)  of  the  chemical 
repeat  unit  constituting  the  polymer  molecule.  However,  there  is 
really  no  theoretical  basis  for  such  a  choice  and,  moreover,  the 
repeat  unit  itself  is  ill-defined  for  some  polymers  (for  example, 
polyethylene).  This  introduces  a  degree  of  arbitrariness  into 
the  numerical  results  calculated  from  the  theory.  The  solution 
to  this  dilemma  is  found  by  dispensing  with  the  use  of  the 
concept  of  segment  size  altogether  and,  instead,  by  recasting  the 
theoretical  expressions  in  terms  of  quantities  which  can 
experimentally  be  determined  unambiguously. 

We  propose  that  the  configurational  properties  of  polymer 
chains  be  expressed,  not  in  terms  of  N  and  a,  but  instead  in 
terms  of  the  molecular  volume  v  and  the  rms  end-to-end  distance  r 


t/- 
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of  the  polymer  molecules.  These  two  sets  of  parameters  can  be 
related  to  each  other  (for  gaussian  chains)  by: 
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r2  =  Na2 


v  =  ca- 


^  m *  .  • 
*  •.  .  •*  .  < 


where  c  is  a  geometrical  constant,  usually  taken  equal  to  1. 

The  Flory-Huggins  free  energy  of  mixing  can  then  be  written, 
in  the  same  spirit,  as  follows: 


AGm  s  kT[ ( 4> i / v i  )ln<J>i  +  2/ x  2) \r\$  2]  +  A(f>^<$>2  (3) 

where  AGm  is  per  unit  volume  of  the  mixture,  4>  -j , 4> 2  are  the  volume 
fractions,  vf,v2  are  the  molecular  volumes,  and  A  is  the 
interaction  energy  density  between  components  1  and  2.  The 
advantage  of  eq.  (3)  over  the  more  traditional  expression  was 
pointed  out  previously . ^ ®  Note  that  the  numbers  of  segments  N ^ 
and  N2  in  the  component  molecules  no  longer  appear  in  (3).  The 
interaction  energy  density  A  is  related  to  the  usual  X  parameter 
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X  =  AVrfif/kT  U) 

where  the  volume  of  reference  Vrej.  is  customarily  taken  equal  to 
the  volume  of  a  solvent  molecule  (if  one  is  involved)  or  some 
average  of  volumes  of  repeat  units  of  the  two  polymers  involved. 
Obviously,  numerical  values  of  x  are  meaningful  only  if  the 
reference  volume  is  clearly  defined  and  specified. 

We  now  rewrite  some  of  the  theoretical  expressions  given  by 
LOW, ^  in  order  to  display  the  explicit  expressions  utilized  for 
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the  numerical  computations  in  this  work.  In  so  doing,  we  also 
incorporate  the  consequence  of  minor  modifications  to  the 
physical  model  originally  employed.  Whereas  LOW^  assumed  that 
homopolymer  A  is  not  penetrating  into  the  micelle  core  consisting 
of  the  B  blocks  of  the  copolymer,  we  now  allow  swelling  of  the 
core  by  the  homopolymer  A  and  seek  to  determine  the  equilibrium 
degree  of  core  swelling  dictated  by  the  minimization  of  free 
energy.  LOW^  presented  explicit  expressions  only  for  the  case  of 
symmetric  diblock  copolymers,  whereas  we  now  write  expressions 
applicable  also  to  non-symmetric  copolymers  in  order  to  be  able 
to  make  comparison  with  experimental  results  given  in  Parts  I  and 
II. 

The  system  of  study  consists  of  a  mixture  containing 
homopolymer  A  of  molecular  volume  vh  and  diblock  copolymer  A-B  of 
overall  molecular  volume  v.  The  volume  fraction  of  A  (or  B) 
units  in  the  copolymer  is  equal  to  fA  (or  fg)  so  that  the  volume 
of  an  A  block  is  vA  =  fAv  (vB  =  fgV,  f A  +  fg  =  1 ) .  We  consider  a 
spherical  micelle  of  radius  R  containing  p  copolymer  chains.  The 
core,  of  radius  Rgf  contains  B  blocks  of  the  copolymer  and  some 
number  of  homopolymer  chains,  so  that  the  volume  fraction  of  B 
monomers  in  it  is  equal  to  nB.  The  spherical  outer  shell  or 
corona,  of  thickness  RA  =  R  -  Rg,  contains  A  blocks  of  the 
copolymer  and  some  homopolymer  chains.  The  corona  contains  only 
A  monomers,  and  the  volume  fraction  of  these  monomers  that  belong 
to  copolymers  is  equal  to  nA. 

The  free  energy  of  a  single  micelle  consists  of  the 
following  contributions: 


(5) 
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where  y  is  the  interfacial  tension  at  the  boundary  between  the 
core  and  corona,  Fd  is  the  contribution  due  to  deformation  of 
copolylmer  chains,  and  FmA  and  FmB  are  the  free  energies  of 
mixing  the  homopolymer  molecules  with  the  blocks  of  copolymer  in 
the  corona  and  the  core,  respectively.  The  interfacial  tension 
can  be  written: 

Y  =  Y0nB  <6) 

where  yB  is  the  interfacial  tension  between  the  corona  and  the 
unswollen  core  containing  only  the  B  blocks  of  the  copolymer. 
Assuming  that  Yq  is  equal  to  the  interfacial  tension  between  two 
homopolymers  A  and  B,  we  can  write  (see  Appendix): 


-Q  -  i/A_\ 

kT  '  3  \kT/ 


1/2 


(r.2/v.)  +  (r*2/v.)1/2(rn2/vB)1/2  +  (r  2/vn) 


(rA"/V 


_ B  '  B‘ 

1/2  ♦  (rB2/vB)1/2 


(7) 


The  contributions  Frf ,  FmA,  and  FmB  can  be  written  as: 

Fd  =  (3/2)  kT  p  (Ra2/«"a2  +  rA2/RA2  +  RB2/rB2  +  rB2/RB2  “  4  ^  (8) 

FmA  =  (4/3)tt  (R3-Rb3)  (kT/vh)  (l-nA)  ln(l-nA)  (9) 

FmB  =  (^/3)^Rb3  t(kT/vh)(i-nB)in(1-nB)  +  AriB(  1  —  nB ) ]  (10) 


From  the  incompressibility  conditions  we  also  have: 


( M/3 )tt  RB3nB  =  pvB 


(11) 


( 4/3 ) tt ( R3  -  R  3)n 


PV  A 


(12) 


Combining  eqs.  (5)— (12),  we  have: 


F/kT  =  pf  =  B.,(  Y/kT)(pvB/nB)2/3 

+  B2(l+c)p5/3(vB/nB)2/3/rB2 

+  B3(1+1/C)p1/3( nB/vB)2/3rB2 

+  P  {( vA/vh)[  <  1-hA)/nA]  ln(1-nA) 

+  (vB/vhH  (  ln(l-ng) 

+  ( A/kT)vB(l-nB)  -  6}  (13) 

where  B1  =  tt1/362/3,  B2  =  35/3/27/3*2/3,  ^  =  61/3*2/3,  and 

C  =  (rB2/rA2)[(vAnB/vB  V1)1/3  "  1]2  (U) 

Next  we  consider  the  free  energy  FM  of  the  whole  system  of 
volume  V.  The  overall  volume  fraction  of  the  copolymer  is  <J>,  and 
of  these  only  fraction  £  participates  in  the  micelles,  the 
remainder  being  molecularly  dissolved  in  the  continuous 
homopolymer  phase.  We  then  write: 

fh  =  CV*?/v>m  ♦  Fmlx  -  TSm  (15) 

where  Fmix  is  the  free  energy  of  mixing  homopolymers  and 
copolymers  in  the  bulk  phase  outside  the  micelles  and  Sm  is  the 
translational  entropy  of  the  gas  of  micelles.  They  are  given  by: 

Fmix/kT  =  V(1-£K>  ((♦1/v)ln*1 
+  C  ( 1- 4>i  J/v^IlnC  1- ) 

+  (A/kTHjfgd-^fg)}  (16) 

and 

Sm/k  =  -V  { (<t^/pv)ln(C4>C  ) 

+  [(1-£<K)/Upv)]ln(1-€<K)}  (17) 


where 


C  =  (1/v)(vA/nA+vB/nB) 


(18) 


and 

♦1  =  <f>(  1-C)/(  1-SH)  .  (19) 

The  equilibrium  values  of  the  parameters  p,  nA,  nB,  and  C  can 
then  be  obtained  by  minimizing  with  respect  to  these 
quantities  with  <J>  fixed. 

NUMERICAL  RESULTS  AND  DISCUSSION 

The  value  of  A  (in  cal/g)  as  a  function  of  temperature  was 
calculated  by: 

A  =  0.71 8±0 .05 1  -  (0.0021±0.000H5)(T°C-150°C)  (20) 

which  was  obtained  by  averaging  the  values  determined  previously 
from  the  cloud  points  of  mixtures  containing  a  polystyrene  and  a 
poly  butadiene  or  a  styrene-butadiene  copolymer.  (The  error 
limits  denote  the  standard  deviation.)  The  molecular  volumes  of 
the  polymer  chains  were  calculated  from  the  molecular  weight  by 
using  the  following  values  of  specific  volume: 

vA  =  1.1138  +  8 . 24x  10_1*t  (21) 

for  polybutadienes2  and 

v  B  =  0.9217  +  5 .412x10“^t  +  1.687x10"7t2  (22) 

for  polystyrenes,2'7  where  vA  and  vB  are  in  cm®/g  and  t  in  °C. 
The  end-to-end  distances  were  calculated  from  the  molecular 
weight  by  r  =  const  x  M^2  with  the  constant®  equal  to  0.9A  for 


poly  butadiene  and  0.7A  for  polystyrene.  The  number-average 
molecular  weights  of  the  polymers,  utilized  for  the  computation, 
are  listed  in  Table  1.  The  minimization  of  the  free  energy  FM 
with  respect  to  p,  n^,  ng,  and  c  was  performed  on  a  PDP  11/31* 
minicomputer.  (The  FORTRAN  program  for  this  computation  can  be 
made  available  on  request.) 

In  Figure  1  the  calculated  values  of  the  core  radius  <R>  are 
compared  with  the  observed  z-average  radius  <R>z.  The  agreement 
between  the  observed  and  calculated  values  is  very  good.  The 
theory  ignores  the  possible  polydispersi ty  of  the  core  radius  (in 
fact,  it  predicts  the  polydisper sity  to  be  very  small),  so  that 
the  calculated  <R>  does  not  distinguish  different  types  of 
averages.  If  we  had  evaluated  the  number-average  <R>n  instead  of 
the  z-average  <R>z  from  the  experimental  data,  the  agreement 
would  probably  have  been  better.  The  predicted  temperature 
dependency  of  the  core  radius  <R>  is  very  small,  whereas  the 
experimental  data  show  that  over  a  temperature  interval  just 
below  the  dissolution  temperature  the  observed  radius  <R>Z  varies 
appreciably  and  goes  through  a  minimum.  The  theory  was  developed 
to  apply  for  the  cases  where  the  incompatibility  between  A  and  B 
units  is  high,  and  it  is  therefore  not  surprising  that  the 
features  observed  near  the  cissolution  temperature  is  not 
predicted  in  detail. 

In  Figure  2  the  observed  and  calculated  values  of  the  CMC 
are  compared.  In  Figure  3  the  same  information  is  replotted, 
this  time  with  the  concentration  on  a  logarithmic  scale.  The 
trend  of  increasing  CMC  with  increasing  compatibility  between  the 


copolymer  and  homopolymer,  as  the  proportion  of  styrene  block 
decreases,  is  correctly  displayed  by  the  calculated  curves.  The 
theory  also  correctly  reproduces  the  general  trend  of  the  CMC 
increasing  very  rapidly  with  temperature  before  approaching  the 
final  dissolution  temperature.  Quantitative  agreement  is  less 
satisfactory,  and  this  is  more  apparent  in  Figure  3  where  the 
disagreement  at  low  concentrations  is  greatly  exaggerated. 

The  most  important  discrepancy  between  the  observation  and 
the  prediction  is  with  respect  to  the  degree  of  swelling  of 
styrene  cores  by  the  butadiene  homopolymer.  Figure  7  of  Part  I 
and  Figure  6  of  Part  II  show  that  the  volume  fraction  ti  of 
styrene  in  the  core  (n  there  is  equivalent  to  nB  of  the  present 
paper)  decreases  steadily  from  unity  as  the  temperature  is  raised 
above  a  certain  temperature  (which  in  the  case  of  sample  75/25  is 
about  90°C).  The  value  of  nB  predicted  from  the  theory,  however, 
falls  from  unity  to  at  most  0.96  with  increasing  temperature 
throughout  the  temperature  ranges  studied  experimentally.  The 
prediction  of  the  very  low  degree  of  swelling  of  styrene  core  by 
poly  butadiene  basically  reaffirms  the  incompatibility  between  the 
styrene  and  butadiene  units,  which  was  deduced  from  the 
miscibility  study^  of  mixtures  of  polystyrene  and  polybutadiene. 
The  value  of  n  falling  as  low  as  0.3  before  the  micelles  finally 
dissolve,  as  found  experimentally,  cannot  therefore  be  explained 
within  the  framework  of  the  model  on  which  the  LOW  theory  is 
based.  To  explain  the  observed  values  of  n.  some  modifications 
may  have  to  be  introduced  into  the  model,  allowing,  for  example, 
the  possibility  of  a  butadiene-rich  inner  core  developing  within 


the  sytrene-rich  core,  or  the  possibility  of  the  micelle  core 
adopting  shapes  different  from  spheres. 

Despite  some  of  the  discrepancies  noted  above,  it  would  be 
fair  to  say  that  LOW  theory  pr  .aicts  the  semi-quantitative 
behavior  of  the  block  copolymer  micelles  well,  especially  at 
temperatures  far  below  the  dissolution  temperature,  and  should 
serve  adequately  as  a  basis  for  understanding  the  thermodynamics 
of  mixtures  containing  a  block  copolymer  and  a  homopolymer. 

After  this  work  was  completed,  a  paper  by  Whitmore  and 
Noolandi^  was  published,  in  which  they  extend  the  earlier  theory 
of  block  copolymer  micelles  by  Noolandi  and  Hong^  and  compare 
its  predictions  with  experimental  results  obtained  by  Selb 
et  al.11  on  block  copolymer-homopolymer  blends  very  similar  to 
ours.  Their  numerical  results  were  calculated  on  the  basis  of 
the  concept  of  a  segment  equal  to  Kuhn  statistical  length,  and 
therefore  cannot  be  compared  directly  with  our  results.  Their 
calculated  values  of  micelle  core  radii,  however,  agree 
approximately  with  our  calculated  values,  and  the  overall  trend 
of  the  CMC  calculated  by  them  also  agrees  with  the  prediction 
from  LOW  theory.  In  view  of  the  close  similarity  between  the 
models  used  in  these  two  theories,  we  believe  that  the  degrees  of 
agreement  with  our  experimental  results  offered  by  them  are 
likely  to  be  comparable. 
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1.  The  radii  <R>  of  the  micelle  core  calculated  from  the 
theory  (thick  lines)  are  compared  with  the  z-average 
radii  <R>Z  obtained  experimentally:  solid  circles  - 
3%  sample  75/25;  open  squares  -  sample  50/50;  solid 
triangles  -  Ut  sample  25/75.  The  theoretical  values 
were  calculated  for  the  same  concentrations  of  the 
copolymers,  but  the  concentration  dependence  is  small 
in  any  case. 

2.  The  observed  and  calculated  values  of  the  critical 
micelle  concentration,  as  a  function  of  temperature, 
are  given  for  the  three  copolymers. 

3.  The  same  data  as  in  Figure  2  are  replotted  with  the 
concentration  axis  on  a  logarithmic  scale. 


ii» 


APPENDIX 


Here  we  give  a  brief  derivation  of  the  interfacial  tension  Yq 
and  the  thickness  d^pt  of  the  interfacial  layer  between  two 
polymers  in  the  long  chain  limit  and  phrasing  the  expressions  in 
terms  of  v,  r  and  A  rather  than  the  customary  set  N,  a  and  x* 

Near  the  interface  the  local  free  energy  density  (per  unit 
volume)  can  be  given  by: 

F/kT  =  (4>/v1  )in  4>  +  t<  1-«J>)/v23ln(1-d>) 

+  (A /kT)  $(  1  -  d>)  +  gradient  term  (A1) 

To  derive  the  gradient  term,  we  consider  the  result  of  random  phase 
approximation1 2  applied  to  polymer  mixtures: 

(s(q)]'1  =  [<t>  vn  fD(  x  n )  3”1  +  EC  1-d>)v2fD(x2>]“1  -  2A/kT  (A2) 

where  sT(q)  is  the  Fourier  transform  of  the  correlation  function,  fD 
is  the  Debye  function: 

fD(x)  =  2(e“x+x-D/x2  (A3) 

and 

x1  =  q2r12/6  (AH) 

For  x>>1 ,  fD(x)  =  12/q2r2,  so  that: 

[s(q)]_1  a  ( q2/1 2 )  [  (  r^/v^  )/0  +  (  r22/ v2)/ ( 1  -  <}>)  ]  -  2  A/kT  (A5) 

The  gradient  term  in  (A1)  comes  from  ( 1/2)^[s(q)]_1  '  6 4>( q )  6 C  — q ) » 
and  utilizing  (A5)  we  find: 

gradient  term  =  (1/24)  [(r.,2/v. |)/4>  +  (»*22/v2)/(  1  —  <J> )  ]  ( Vd> ) 2  (A6) 


In  the  limit  of  long  chains,  (A1)  can  therefore  be  written  as: 

F/kT  =  (A/kTH  ( 1-<|> )  +  ( 1/24 )  [  ( r  1 2/ v1  )/4>  +  (  r22/ v2 )/ ( 1- <p)  ]  (V<*>  )2  (A7) 

Next  we  derive  the  density  profile  at  the  interface  boundary. 
Applying  the  Lagrange-Euler  relation: 

9  F/34>  =  (9/9  x ) ( 9  F/3$ )  (A8] 

(where  4>  =  9$/9x)  to  (A7),  we  find: 

(A/kT) ( 1-2  <p)  =  (<P/12)[(r12/v1)/^  +  (r22/v2)/n-<t,)] 

+  (&2/2 4)[-(rj2/v.j  )/4>2+(r22/v2)/(  1- $>)2]  .  (A9‘ 

Multiplying  both  sides  of  (A9)  by  4>  and  integrating  with  respect  to 
x,  we  obtain  (for  a  very  incompatible  case  where  4>-»- 0  and<J>->-0  or  1  as 


x  -*•  ±oo)  : 


(A/kT)4>(1-$)  =  (4>2/2i J)  [  ( r  1  v1  )/<J>  +  ( r2c/^2)/ ( 1-4>)  3 


\-t4 


(A10) 


The  density  profile  can  be  derived  by  integrating  34>/3x  obtainable 
from  (A10).  In  particular,  the  thickness  dint  of  the  interface 
layer,  defined  as  the  reciprocal  of  the  gradient  at  41  =  1/2,  is  found 
to  be: 


r,2/v,  .  r22/v; 


(AID 


The  interfacial  tension  yQ  is  obtained  by  integrating  the  free  energy 
F  given  by  (A7)  (with  the  density  profile  given  by  A10  substituted  in 
it)  for  the  whole  range  of  x  between  -  «  and  °°  .  The  result  is: 

Y0/kT  =  (A/kT)1/2(2/3)[(r12/v1)3/2-(r22/v2)3/2]/(r12/v1-r22/v2).  (A12) 


'  •  *  *  '  *  . . 
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